%**************************************************************************
%
%                           估价函数
%
%**************************************************************************
sr = [Sr(1, K+1: length(Sr)), zeros(1,K)];
delta = st - sr;
Perr = sum(delta.^2) / N;
U = K + 10 * log10(Perr);

